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Abstract. We present a new matrix- valued isospectral ordinary differen- 
tial equation that asymptotically block-diagonalizes a finite-dimensional zero- 
diagonal Jacobi matrix employed as its initial condition. This differential equa- 
tion is closely related to the one introduced by M. Kac and P. Van Moerbeke 
in 1975, although our approach to prove the key properties of this o.d.e. dif- 
fers from the techniques developed by them. We show that our o.d.e. can be 
represented as a double bracket differential equation similar to the one studied 
by R.W. Brockett in 1991. 



§ 1. Introduction and Main Result 

Solutions to tasks of sorting a list, diagonalizing matrices, and solving linear pro- 
gramming problems are traditionally associated with computer science algorithms, 
for example the quicksort algorithm for sorting or the simplex method for solving 
linear programs. Brockett [Brockett, 1991] showed that solutions to such problems 
can also be obtained by means of a smooth dynamical system, in particular, as the 
limit of solutions to certain matrix- valued ordinary differential equations (o.d.e.'s). 
A classical problem from linear algebra is therefore solvable by calculus. Motivated 
by Brockett's work, new problems, conventionally tackled by algebraic methods, 
have been assigned to calculus. For instance, [Faybusovich, 1992] proposed an or- 
dinary differential equation (structurally similar to the one proposed by Brockett) 
as the starting point in a general approach to interior point methods for linear 
programming. 

By a Jacobi matrix we mean a symmetric tridiagonal matrix (in general, infinite) 
with real entries and distinct eigenvalues. In this article we present a matrix- 
valued ordinary differential equation which asymptotically block-diagonalizes a 
finite-dimensional zero-diagonal Jacobi matrix taken as its initial condition. Ja- 
cobi matrices arise in a variety of applications, for example in solid state physics 
to characterize the Toda lattice, which is a simple model for a one dimensional 
crystal — see e.g. [Moser, 1975] for a detailed study. Indeed, the Toda lattice is a 
dynamical system consisting of n particles in a potential field described in terms of 
the Hamiltonian 

H(qi, . . . , qn,pi, ■ . ■ ,Pn) '■— Potential Energy + Kinetic Energy 
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n— 1 n 

m— 1 m— 1 

The particle m, therefore, interacts only with particle m — 1 and particle m + 1, 
but the positions of the particles are not required to be ordered. Also observe that 
the potential energy X]m=i e''"~'™+'^ is asymmetric, and for a given level of energy 
E (kinetic + potential), one has — Qm+i ^ InE. In an appropriately chosen 
coordinate system, the dynamics of this completely characterized by a matrix o.d.e. 
involving a Jacobi matrix; for details see, e.g., [Hclmke and Moore, 1994, pp. 59-60]. 
There is also a strong connection between Brockett's double bracket flow [Brockett, 
1991] and the Toda lattice equation, which was first observed by [Bloch, 1990]. 

We offer a second motivation here that has intrinsic appeal and interest, relating 
to the computation of the roots of certain polynomials. Orthogonal polynomials 
on the real line corresponding to a Borel probability measure /i have considerable 
applications in mathematical physics and engineering [Simon, 2005]. Let (•,•) de- 
note the standard inner product on the Hilbert space L2(M,/^). Then the monic 
orthogonal polynomials on the real line are defined recursively [Szego, 1959] by 

Pn+l{x) = xP„{x) - a^,P„_i(x) - bn+lPnix), U G N, 

P_i(-) = 0, Po(-) = l, 

where 

ll^'^ll (xP.n,Pn) 

a„ := 177- — \, for n e N and b^+i ^ ^ for n £ Nq. 

Il^n-lll \\Pn\\ 

To a given measure /z, there corresponds a (generally infinite) Jacobi matrix 

/bi ai •••\ 

fli 62 0.2 • • • 

J 02 63 03 • • • 

03 &4 ••• 
V: : : ; 

with strict positive off-diagonal entries derived from its orthogonal polynomials, and 
it is well-known that the zeros of these orthogonal polynomials are precisely the 
eigenvalues of finite truncations of this Jacobi matrix. Conversely, Favard's The- 
orem [Favard, 1935] shows that to every finite-dimensional symmetric tridiagonal 
matrix with strictly positive off-diagonal entries there corresponds a finitely sup- 
ported measure. (The uniqueness of this measure for an infinite-dimensional Jacobi 
matrix is an issue that relates to solvability of the "moment problem" [Akhiezer, 
1965], with which we shall not deal here.) Since the space of square integrable func- 
tions corresponding to this measure is finite dimensional (the measure itself being 
finitely supported) , it is enough that from the orthogonal polynomials it is possible 
to recover the measure, and this leads to the problem of finding the roots of such 
polynomials. The latter, in general, is known to be a difficult task. As mentioned 
above, for our applications it suffices to determine the eigenvalues of the finite trun- 
cations of the corresponding Jacobi matrix, which we shall do in this article with 
the aid of an appropriate matrix-valued ordinary differential equation. For this we 
further specialize the measures to non-negative linear combinations of finitely many 
Dirac measures on the real line placed symmetrically around 0. These measures 
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give rise to zero-diagonal Jacobi matrices, and our main result to the problem of 
finding the roots of the orthogonal polynomials corresponding to these measures. 

To summarize, the problem that we treat in this article involves is that of com- 
puting the eigenvalues of zero-diagonal, finite Jacobi matrices. 

Some preliminary notation follows. We let Sym(n) and Skew(7i) be respectively 
the set of symmetric and skew-symmetric n x n matrices with real entries. We 
define Jaco(fi) as the set of all n x n Jacobi matrices with real entries and zeros 
on its diagonal. For a matrix A, \\A\\ is the Frobenius norm defined as ||j4|| = 
^yti^AA^). The bracket [•, •] is the usual matrix commutator [A, B] = AB — BA. 
For Hq e Sym(n) we let M(i?o) {O^ HoO \ 6 e 0{n,R)} denote the set of all 
real matrices orthogonally similar to Hq [Moore et al., 1994]. For Hq e Jaco(n) we 
define Mj(_ffo) to be the set of all zero-diagonal Jacobi matrices that are isospectral 
to Hq (that is, they have the same eigenvalues). Moreover, for p — 0, 1, . . . , n — 1, 
we let 



i-p) 











ai 




02 



^n—p 





V 



/ 

Vp{ai,a2, . . . , CLn-p) '■— T^u,p{0'l, 02, ■ • • , On-p) + ^^u,p{a,l^ CL2: ■ ■ ■ : a„_p)^. 

(1.1). Theorem. Let n > 1 be given. Consider the zero-diagonal Jacobi matrix 

\ 



/O ai 

fli 02 

0.2 



(1.2) 



On- 2 




On-2 

a„_i 
o„-i 



e Jaco(n), 



(1.3) 



and the skew-symmetric matrix 

K{H) :=I?„,2(ai02,02a3, ■ ■ • ,a„_2a„-.i) 

— I'n, 2(0102, O2O3, • • • , a„_20„_i)^, 

derived from H. Consider the matrix-valued o.d.e. 
d 



(1.4) 



dt 



H{t) = [H{t),K{H{t))], 



H{Q) := Hq e Jaco(n). 



(i) (1.4) defines an isospectral flow on the set of all Jacobi matrices with zero 
diagonal entries. 

(ii) The solutions H(t) of (1.4) exist for allt^ and approach asymptotically the 
set of equilibrium points of (1.4), which is given by 

{H eMj{Ho)\K{H) = 0}. 
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(iii) limi_5,oo H{t) exists and it is a block diagonal matrix, where each block is of 
size 2x2. If H(t) has an odd size, there is exactly one 1x1 block containing 
a zero. 



1. Brockett [Brockctt, 1991] studied the o.d.e. (1.4) in which the map K on the 
right-hand side of (1.4) was set to K{H) = [H,N], where is a constant 
symmetric matrix. On the one hand, we shall see that in our case the role of 
N is played by a linear function of H, not a constant; as such, the analysis of 
(1.4) requires new tools. On the other hand, whereas it is possible to treat N 
as a parameter in [Brockett, 1991], and can therefore be employed, e.g., to sort 
the eigenvalues of limt_>.oo H(t) according to a particular order by selecting an 
appropriate iV, such variations in our case are not readily available. 

2. Although the dynamical system -^i?(t) = [K {H {t)) , H {t)] was studied first 
in [Kac and van Moerbeke, 1975], their techniques were entirely different. In 
particular, the analysis of its properties from a double bracket perspective, to 
our knowledge, has been carried out here for the first time. 

Outline of the article: §2 presents the proof of the Theorem, which we illustrate 
with some numerical examples in §3. We conclude in §4 with presenting an extension 
of our work and discuss possible channels of further research. 



Some preliminaries are needed in order to prove Theorem (1.1). We begin with 
the following classical result, which will play a key role behind proving that the 
solutions to (1.4) are isospectral. 



Then there exists a smooth family of unitary matrices (C/(i))t^o with U(0) — /, 
such that 



Remarks. 



§2. Proof of Theorem (1.1) 




Skew(n) be a smooth map 



j{t) = u{ty'^{o)u{t), 

The family (7(i))t^o is thus isospectral. 



t ^ 0. 



Next we define the mapping 



(2.3) 



/an ai2 



\ 



Sym(n) 3 A 



'^n— l,n— 1 ^n — l,n 



N{A) := X'i(-ai2, 0, 034, 2045, . . . , (n - 3)a„_i 



,n 



) e Sym(n). 
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The mapping N(-) is linear, and can be written as 



n-l 



(2.4) 



where A e Sym(n) and Ei is the matrix with 1 at its (i, i)-th entry and zeros 
elsewhere. 

(2.5). Proposition. Let 



/O ai 

fli a2 

02 



(2.6) i/„ 



V 



a«-2 
a„_i / 



e Jaco(?T.). 



Then the commutator of H„ and N(Hn) is given by 

^ ^ K{H^):=[H„,N{Hn)] 
(2.7) 

= I'u, 2(01^2, • ■ ■ , an-2an-l) — 2'«,2(aia2, • ■ • , O'n-20'n-l) , 

where N{-) is the linear mapping defined in (2.3). 



Proof. The proof proceeds by induction. Observe that I?„ 2 = if n = 1. For n = 1, 
we have Hn = 0, and therefore K{Hn) = 0. For n = 2, we have _ff„ = ( ) and 
N{Hn) = („°^ ). Since H,, and 7V(i/„) commute, X(i7„) = 0. 

We consider the induction step: 



(2. 



Hn+l — 










o\ 




Hn-l 























... 





a-n-i 


a„ 


Vo ... 





a„ 


0/ 



and 
(2.9) 





VO 



iV(i7„ 



[n - 3)a„_i 








(n — 3)a„_i 

(n - 2)a„ 

(n - 2)a„ y 
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Simple matrix multiplications lead to 
(2.10) 



Hn+lN {Hn+l) 



where 



and 
(2.11) 



/ 

















i/„_iiV(i/„_i) 



















u 








(n 


— 3)a„_ia„_2 
















® 





... 


(n - 4)a„_ia„_2 





® 





\o 


... 





® 





@/ 




® (n - 3)a„a„_i, 










®=in-2)al 










(3) = (n - 2)a„a„_i, 












-2)a^, 






/ 











o\ 






7V(i/„_i)i/„_i 




























(n 


- 4)a„_ia„_2 
















® 





... 


(n - 3)a„_ia„_2 





® 





Vo 


... 





® 









® 


= (n - 2)a„a„_ 










® 


= (n - 3)a„a„_ 


1- 







where 



By the induction hypothesis, 

[H„^i,N{H„^i)] = i/„_iiV(i7„„i) - N{Hn-i)Hr,-i - K{H„^i). 
Therefore, by (2.10) and (2.11) we get 



[F„+i,iV(iJ„+i)] 











\ 





an-ian-2 

a„a„_i 

... -a„_ia„_2 

Vo ... -a„a„_i / 

if(if„+i). 



□ 
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Note, that in view of Proposition (2.5), the modified Kac-van Moerbeke equation 
(1.4) can be represented as the double bracket o.d.e.: 

(2.12) -^^Hit)^[Hit),[H{t),N{Hm]. i^o := i^(O) G Jaco(n). 
We shall employ the following auxiliary lemma in the proof of Theorem (1.1). 

(2.13) . Lemma. // A = I?i(ai, a2, a„_i) with ai R for all i, then 

[A,[A,N{A)]] = 

\ / ^ /222 22 2^ 

T>i[—aia2, —0203 + a^a2, ■ ■ ■ , — a„_2a„_i + a„_3a„_2, a„-2'^n-i); 
where the mapping N(-) is as defined in (2.3). 

Proof We have to show that K{A) = [A,N{A)]. To this end, let Ai denote the 
upper triangular part and A2 the lower triangular part of ^, such that A — A1 + A2. 
Using an analogous decomposition for K{A) we get Ki{A) + K2{A) = K{A). We 
abbreviate and simply write Ki for Ki{A), i G {1,2}. We decompose 

[A,K{A)] = [Ai + A2,Ki+K2] 

= [AuKi] + [Ai,K2] + [A2,Ki] + [A2,K2]. 

We first show by induction that [Ai,iiri] = 0. Indeed, if we denote with a su- 
perscript the size of a matrix and consider = I?i(ai, 02, . . . , On-i), i-e.. A" ~ 
2^«,i(ai, «2, • • • jfln-i) and A2 = X'„_i(ai,a2, . . . ,a„_i)^, we observe that for n = 1 
we have [A", K^] — 0. The induction step can be done as follows: 



( 



Vo 



A'l 



\ 





an 

/ Vo 

\ / 



(^n— 1 f^n 



/ 



dfi— 1 



A'^ 



I Vo 
\ 



Vo 















an 
/ 



Vo 














= 0. 



Moreover, we have A^ ~ A2 and KJ = — i^2- Therefore, 
[A2, K2] = A2K2 - K2A2 = -AlKl + KJAI = {A^K^ - K,A,y = [A,J<,Y 
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Thus, [A2, K2] = as well, and it remains to show (again by induction) that 

[A2,Ki] + [A^,K2] 



= 'Di{—aia\, —0203 + a\a2 



, . . . , U/j^. 



2,2 2 N 

-2a„_i + a„_3a„_2, a„_2a„_ij. 



As the next step, we claim that 

[^2 , X"] = 'Du,i{-aia\, -a2al + a\a2, . • • , -a„_2a^-i + a^_3a„_2, a?j-2«"-i)- 
The induction base is trivial. Then we have: 



AS 






o„ 0/ 



\0 ... 

\ / 



\0 



Vo 





/ 
\ 





al_ia, 
/ 



Vo 





/ 




a„ 0/ 



^2 







fln-ia^i 



VO ... 0/ yo ... 0/ 

= T^u^i-a-ial, -0203 + a^a2, . . . , -a„_ia^ + a^_2an-i, a^-iOn)- 
Again, in view of = yl2 and KJ — ~K2, we obtain 

[Ai,K2] = A1K2 - K2A1 = -AlKl + KJAI = {A2K1 - KiA2Y = [A2,KiY . 
Thus, 

[A", i^a'] = 'Du,i{aial, a2al + a\a2, . . . , -a„_2a^_i + a^_3a„_2, a^-2an-i)^- 
According to (2.15), for A = I?i(ai,a2, . . . ,a„_i) we have 

[A, K{A)\ = 'Di{-aia\,-a2a\ + 0^02, . . . , -a„_2a^_i + a\_^an-2,a'^-2<^n-i), 
which completes the proof. □ 
(2.16). Lemma. For A,B £ Sym(n), 

[A, [A, B]\ = if and only if [A, B] = 0. 
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Proof. The "if" part is trivial. To prove the "only if" part, suppose that [A, [A, B]] — 
0. This implies B[A, [A, B]] = 0. Using the techniques in [Helmke and Moore, 1994, 
p. 49], we compute 

0-tr(B[A, [A,B]]) 
= tr {B{A^B - 2ABA + BA^)) 
= tr {BA^B - 2BABA + B^A^) 
= tr {BA^B - BAB A - ABAB + AB'^A) 
= tr ((BA - AB){AB - BA)) 
= tr([B,A][A, B]) 
= tr([A,B]-[A,B]) 

which immediately gives [A, _B] = 0. □ 

(2.17) . Lemma. Consider the continuous Junction 

(2.18) M{Ho)3H^f{H) := -^\\H - NiH)f +^-\\NiH)f eR. 
With respect to the o.d.e. 

(2.19) ^^Hit) = [H{t),[H{t),N{Hm], H (0) ^ Ho e Sym{n) , 
the time derivative of f{H{-)) is given by 

(2.20) ^^f{H{t))^\\[H{t),N{H{t))W. 

Proof. We start by simplifying the function 
f{H{t)) = -\\\H{t) - N{H{t))r + \\\N{Hit))r 

= tr {H{t)H{t) - H{t)N{H{t)) - N{H{t))H{t) + N{H{t))N{H{t))) 
+ \\\N{H{t))r 

= -\mt)r + \tr{N{H{t))H{t)) - \\\N{H{t))r + \\\N{H{t))r 

= -\\\H{t)r + \tr{N{H{t))H{t)). 

Since Hit) e M(i?o) by Proposition (2.1), \\H{t)\\ is constant for all < ^ 0. We 
calculate the derivative of / along trajectories of (2.19) as follows: 

^J{H{t)) = \tr {J^±N{H{t))^ H{t)+N{H{t))H{t)^ 

-\^^{{i^im))Hit)) 

+ \tr{N{H{t))[H{t),[H{t),N{Hm]) 



\i.[[lNmt)))Hit)) 
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1 



[N{H{t)),HmH{t),N{Hm 



(2.21) = [(y^N{B{t))\^ H{t)^ + ^-\\[H (t) , N {H , 

where, at the third equahty, we employed the fact [Bernsteni, 2009, p. 162] that for 
A,B,C e M"""", tr{A[B, C]) = tr{[A, B]C). Therefore, it remains to show that 



tr (^(^^iV(i/(i))) H{t)'^ = \\[H{t),N{Hm\\ 



2 



Note that 

-^^N{H{t))^N{H{t)) 

since iV is hnear, and since from (2.21) it follows that tr(7V(i/(t))ij(i)) = \\[H (t) , N {H {t))]\\' 
our proof will be complete if we show that 

(2.22) tr{N{H{t))H) = tr {N {H {t)) H (t)) . 

To this end, employing the expansion of N in (2.4), we see that 



/n-l 

tr {N{H){t)H{t)) = tr ^(z - 2) {E,H{t)E,+, + 

\i=l 


,iH{t)E,) H{t)^ 


= ^(z - 2) tr {{E,H{t)E,+i + E,+ 

i=l 


iH{t)Ei) H{t)) 


= ^(z - 2) (tr {E,H{t)E,+,H{t)) 

1=1 


+ tr{E,+iH{t)E,H{t))) 


= ^(i - 2) (tr {E,+iH{t)E,H{t)) 

1=1 


+ tr {E,H{t)E,+^H{t))) 


ri-l 

= ^(z - 2) tr {{E,H{t)E,+, + E,+^H{t)E,) H{t)) 

A — 1 


^tr{N{H{t))H{t)), 





which establishes (2.22), and completes the proof. □ 

(2.23) . Lemma. The o.d.e. 

(2.24) ^^H{t) = [H{t), [H{t),N{H{t))]] , H{Q) = Ho e Jaco(n), 

has a finite number of equilibrium points in Jaco(n.) that are isospectral to Hq. 

Proof. In view of Proposition (2.1) and Lemma (2.13), H{t) G Mj{Hq) for all 
solutions of (2.24) and for all t ^ 0. Moreover, for Hq S Jaco(7T.) we have Mj(i?o) Q 
Jaco(n). By Lemma (2.16), the set 

E:^{H{t)eM,{Ho)\ \\[H (t) , N {H mf = O} 
= {H{t) e Mj(i7o) I [H{t),N{H{t))] = 0} 
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is the set of all equilibrium points of (2.24) on Jaco(n) that are isospectral to Hq. 
Let 

H = 'Di{ai,a2,a3,ai, . . . ,a„_2,a„_i) S Jaco(n). 

At this point it is crucial to recall that according to our definition Jacobi matrices 
have distinct eigenvalues. We treat the case of n even and n odd separately: 

n even ] 

Consider the set of matrices 

E:^{H e Jaco(n) | [H,N{H)] = O}. 

In view of (2.7), the only possibility for H to lie in is if = for all i even and 
Ui =/= 0, for all i G {1, 3, ... n — 1} with a,; ^ aj for all i =/= j, such that 

He = 2^1(01,0, 03, 0,...,0,a„_i). 

Note that Hj^ has the spectrum (j{H^) = {ioi | i G {1, 3, . . . n — 1}} containing 
only distinct eigenvalues. Now E as defined in (2.25) is a subset of E satisfying 
the isospectral conditions; it is the restriction of E to the set of zero-diagonal 
Jacobi matrices isospectral to Hq, i.e., E = Considering all the possible 

permutations of the for iG{l,3,...,n — 1}, the set E contains (^)! equilibrium 
points on Jaco(ri) that are isospectral to Hq. 

n odd 

First of all, since H = T>i{ai, 02, as, 04, ... , a„_2, dn-i) is a (zero-diagonal) Jacobi 
matrix, we need to evoke the fact [Pcnskoi, 2008] that its spectrum has the form 

(2.26) cr(7?) = {0,±Ai,±A2,...,±A^}, where 7^ for alH ^ j. 
Moreover, if we consider 

E:= {H e Jaco(n) | [H,N{H)] = O}, 
in view of (2.7), H must satisfy 

(2.27) flifli+i = for all i = 1, . . . , n — 2 

in order to lie in E. Wc claim that ai and a„_i are nonzero. Indeed, if we assume 
for the sake of a contradiction that ai = 0, then according to (2.27) a„_i = 0. This, 
however, leads to being an eigenvalue of multiplicity two, therefore contradicting 
(2.26), and vice versa. Hence, a2 = 0, a„_2 = and H contains the diagonal blocks 

( ai '^0 ) ^^'^ ( aH-i "V^ ) ^^^^ ^^'^ ^^^^ diagonal blocks. Furthermore, in 

order to satisfy (2.26) and (2.27) H has to be a block-diagonal matrix containing 
2x2 blocks of the form ^ ^, °q ^ with Uj ^ and one 1x1 block containing a 

zero, where the upper and lower blocks are 2x2 blocks and the block entries are 
distinct (since the eigenvalues of of H need to be distinct). Accordingly, there 
are possibilities to place the 1x1 block in H. As above E = i?|Mj(ffo) 
and considering all possible permutations of the blocks we get that E contains 
equilibrium points on Jaco(n) that are isospectral to Hq. □ 

Proof of Theorem (l-l)- In view of Lemma (2.13) it follows that the right-hand 
side of (1.4) is a symmetric tridiagonal matrix with zero diagonal entries given by 
(2.14). The fact that (1.4) is isospectral is an immediate consequence of Proposition 
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(2.1). Therefore, the flow of (1.4) evolves on the set of zero-diagonal Jacobi matrices 
isospectral to Hq, i.e., H{t) e Mj(i/o) for all t ^ 0. This settles the claim in (i). 

In order to show (ii), note that in view of (i), H{t) e Mj{Ho) for all t ^ 0. 
Since Mj(iJo) is known to be a compact manifold [Tonici, 1984], H{t) exists for all 
t ^ 0. By Lemma (2.16) we see that the set of equilibrium points H of (1.4) is 
given by {H e Mj{Hq) \ [H, N{H)] = O}. To show that H{t) approaches the set of 
equilibrium points, consider the function 

(2.28) Mj(i/o) f{H) := -\\\H{t) - N{Hm? + \\\N{Hm? & K. 
According to Lemma (2.17), 

(2.29) ^^f{H{t)) ^\\[H{t),N{Hmr ^Q. 

We invoke the Krasovskij-LaSalle's Invariance Principle [Khalil, 2002, Theorem 4.4], 
[Vidyasagar, 2002, p. 178]: First, we define the set 



E:=i^H e Mj{Ho) [H, N{H)] = o} 



By Lemma (2.16), E coincides with the set of all equilibrium points of (1.4). There- 
fore, E is an invariant set with respect to (1.4). Second, recall that Mj(ifo) is a 
compact set. Therefore, by the Krasovskij-LaSalle's Invariance Principle, every 
solution {H{t))t^o starting in Mj(_ffo) approaches the set of equilibrium points E 
asymptotically, which proves the claim in (ii). 

To show property (iii) note that we have already shown in part (ii) that every 
solution of (1.4) approaches the set of equilibrium points asymptotically. However, 
according to Lemma (2.23), the number of equilibrium points is finite. By continuity 
of trajectories, therefore, {H{t))t^Q converges to a single equilibrium point, i.e., 
limt_>.oo H(t) exists. Furthermore, we have seen in the proof of Lemma (2.23) that 
each of the equilibrium points has exactly the block diagonal structure stated in 
the theorem, which proves the claim in (iii). □ 



§ 3. Examples 



We illustrate Theorem (1.1) with two numerical examples. To give an idea about 
the rate of convergence of the o.d.e. (1.4), we have solved the o.d.e. numerically and 
plotted the upper diagonal entries of H{t) as functions of time for a chosen intial 
condition. The figures show that the the solution of (1.4) converges rather quickly 
to an equilibrium. 

(3.1) . Example. We start with a 7 x 7 example. Consider an initial condition 

(3.2) iJ(0) = i?o =X'i(10,8,10,5,5,9), 

having the spectrum (t{Hq) = {0, ±7.38, ±10.57, ±15.13}. By solving the o.d.e. 
(1.4) numerically, we see that the transient behavior vanishes well before 1 second 
of simulation. At T = Is we have 

H{T) = X>i(0, 7.38, 0, 10.57, 0, 15.13). 

Figure 1 shows the evolution of the super-diagonal components of H{t) against time 
t. 
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Figure 1. Flow of o.d.e. (1.4) for the initial condition (3.2). 

(3.3) . Example. Consider the following initial condition of dimension 16 

(3.4) H{0) = iJo ==X'i(3,10,15,2,6,ll,5,6,12,5,12,4,ll,3,5) 
with the spectrum 

CT(iJo) = { ± 1-73, ±2.85, ±4.76, ±9.51, 

± 12.22, ±13.59, ±15.83, ±18.20}. 

By solving the o.d.e. (1.4) numerically with this initial condition we see again that 
the transient behavior vanishes well before 1 second of simulation. At T = Is we 
have 

H{T) = I?i(1.73, 0,2.85, 0,4.76, 0,9.51,0, 

12.22,0,13.59,0,15.83,0,18.20). 

Figure 2 shows the evolution of the super-diagonal components of H(t) against time 
t. 

§4. Conclusions and Future Direction 

We have analyzed a modification (1.4) of the Kac-van Moerbeke equation in this 
article. It was demonstrated that this o.d.e. can be represented as a double bracket 
equation, thus establishing a connection to [Brockett, 1991], and we established 
certain key properties of this o.d.e. by employing Krasovskij-LaSalle's Invariance 
Principle as one of the major ingredients. 

The domain of the o.d.e. (1.4) can be expanded to the set of real symmetric ma- 
trices Sym(n). Since M{Ho) for Hq G Sym(n) is again a compact manifold [Hclnike 
and Moore, 1994], assertions (i) and (ii) of Theorem (1.1) hold also for the sym- 
metric case, and the proof proceeds analogously. Extensive simulations lead us to 
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Figure 2. Flow of o.d.e. (1.4) for the initial condition (3.4). 

conjecture that the solutions converge asymptotically to block diagonal matrices, 
as in the case of zero-diagonal Jacobi matrices employed as initial conditions. How- 
ever, a proof for this conjecture is still an open problem — in contrast to the case 
of zero-diagonal Jacobi matrices, there exist infinitely many equilibrium points in 
this case. 

Another direction that appears to be particularly interesting for future reseach is 
the sorting of the diagonal blocks. We know that in Brockett's double bracket flow, 
if the matrix N is chosen in an appropriate way, the diagonal elements in the limit- 
ing solution are sorted [Brockctt, 1991]. The numerical experiments reported in this 
article lead us to conjecture that for almost all initial conditions Hq S Jaco(n) the 
solution {H{t))t^o of the modified Kac-van Moerbeke equation converges asymp- 
totically to a block diagonal Jacobi matrix with zero diagonal entries, where the 
blocks are sorted by increasing absolute value of the eigenvalues. 
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